Wilson Fermions and Axion Electrodynamics in Optical Lattices 
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We show that ultracold Fermi gases in optical superlattices can be used as quantum simulators of different 
types of relativistic lattice fermions in 3+1 dimensions. By exploiting laser-assisted tunneling methods, we find 
the atomic analogue of the so-called naive Dirac fermions, and thus provide a physical realization of the fermion 
doubling problem. Moreover, we show how to implement Wilson fermions, and discuss how their mass can be 
inverted by tuning the laser intensities. In this regime, our atomic gas corresponds to a remarkable phase of 
matter where Maxwell electrodynamics is replaced by axion electrodynamics: a 3D topological insulator. 
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The formulation of relativistic fermions in lattice gauge 
theories (LGTs) [1] is hampered by the fundamental prob- 
lem of species doubling [2], namely, the rise of spurious 
fermions that modify the physics at long wavelengths. To pre- 
vent the abundance of fermion doublers, a suitable tailoring 
of their masses is required, leading to the so-called Wilson 
fermions [3]. A different, but also fundamental, hindrance 
in high-energy physics is the strong CP problem, more pre- 
cisely, the lack of experiments confirming the charge-parity 
violation in strong interactions [4]. To reconcile theory and 
experiment, the Peccei-Quinn mechanism postulates a new 
particle, the axion [5], whose detection still remains elu- 
sive. Interestingly enough, these two seemingly unrelated 
problems turn out to be closely connected. Indeed, Wilson 
fermions with an inverted mass give rise to a certain axion 
background [6, 7]. In this Letter, we suggest to exploit this 
connection to explore axion electrodynamics [8] in a table- 
top experiment of ultracold atoms. The exquisite and gen- 
uine control over quantum systems inherent to ultracold-atom 
technologies [9] allows us to propose the implementation Wil- 
son fermions in optical superlattices (Fig. la). This experi- 
ment would trace a promising route to design, control, and 
probe the rich physics of axions. As a first step, inverting 
the Wilson mass via laser-assisted tunneling enables us to re- 
alize a background axion field 9 = n, which corresponds to 
a new class of unconventional states of matter: 3D topologi- 
cal insulators (TIs) [6, 7, 10, 11]. These gapped phases have 
conducting edges protected by topological order, but respect 
time-reversal symmetry. We show that our proposal consti- 
tutes the first fully-controllable quantum simulator (QS) [12] 
of 3D TIs preserving a general anti-unitary symmetry. Be- 
sides, the space- and time-dependence of the axion field can 
be experimentally tailored to test the magnetoelectric [6], Wit- 
ten [7], or Wormhole effects [13]. We focus instead on a frac- 
tional magnetic capacitor, whose implementation and detec- 
tion are better suited for optical -lattice techniques. Along this 
route, several intermediate, but interesting, phenomena can be 
observed in our QS: i) 3+1 massless Dirac fermions; ii) 3+1 
massive Dirac fermions; and Hi) Wilson fermions. 



We consider a K Fermi gas in an optical superlattice [9], 
focusing in the Zeeman sublevels (i.e. spins) of the F = j 
hyperfine manifold. Below, we show that laser-assisted tun- 
neling methods lead to the effective Hamiltonian (h = 1) 

#eff = Y^L^ C l+vA U rv]x'xC r x + H.C, (1) 
rv TT ' 

where cJ T (c rT ) creates (annihilates) a fermion with spin T at 
site r = mx + ny + lz with m,n,l e {1...L}, and t v ~ 0.1-lkHz 
is the tunneling strength. Here, U rv are operators dressing the 
hopping from r — » r + V, V G {x,y,z}, and we use gaussian 
units. These operators usually rely on spin-dependent optical 
lattices [14]. We use instead spin-independent bichromatic su- 
perlattices (Fig. la), which trap all levels from the F = { J, |} 
manifolds, and allow for lifetimes Ti ~ Is. The optical poten- 
tial V(r) = Volv[cos 2 (7rr v ) + cos 2 (27i;r v )], where V - 50- 
150kHz and we have set the lattice spacing to 1. This yields 
a cubic superlattice with atoms trapped in the minima at zero 
energy (i.e. sites), and secondary minima at AE S \ ~ 50-100 
kHz (i.e. links). The hopping U rV between F = 5 atoms in 
neighboring sites is mediated by a Raman transition to a F = 1 
"bus" level in the intermediate link (Fig. la). 

Let us consider two states |cf,Q, |oc', £'), where a, a' label 
the Zeeman sublevels of the F = {5,5} manifolds, whereas 
label the band index and center-of-mass coordinates 
spanning the lattice sites and links. A two-photon process (a), 
and p, are the frequency and momentum of the 2-th photon), 
after eliminating an excited level, couples these states 

4=££ & a aC c \> a > c Z<* + Hx -> = SC'C^W-™ 

(2) 

where c^ a (c^ a ) creates (annihilates) a fermion in \a, Q, S^i^ 
is the overlap between Wannier wavefunctions in presence of 
a momentum transfer pi — P2: Sqiq = (£'|e~'( p2 ~ Pl )' x |£), and 
co = <x>i — ©2. The remaining part of the coupling depends on 
the light polarization and the atomic internal structure. Thus, 
the formula factors out the contribution of the center-of-mass 
wavefunction and of the internal degrees of freedom. If the 
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transferred momentum is large, we increase the overlap factor 
between neighboring sites and links, and thus the hop- 
ping. Besides, if lasers are far-detuned from this transition, we 
can adiabatically eliminate the "bus" level, realizing a four- 
photon coupling that leads to the single-band Hamiltonian in 
Eq. (1). Here, the hopping strengths scale as t v ~ |£2| 2 /of, 
where d ~ 0.2 - 2MHz, and i2 - 10-100kHz is the Rabi fre- 
quency (Fig. la). Due to the large Zeeman shift AE Z /B ~ 0.3 
MHz/G, we can independently implement each matrix ele- 
ment of U rv eliminating a different "bus" level. A careful 
analysis shows that the contributions of other bands and spu- 
rious on-site couplings can be neglected [15]. In Figs, lc-d, 
we confirm this for two schemes with the necessary ingredi- 
ents: spin-preserving and spin-flipping hoppings. We use 4 
spin components, and design the hopping in terms of Pauli 
matrices (Fig. lb), U rV — e _1( ^ vav , where a v — a z ® O v and 
(p v € M. Such block-structure allows the implementation of 
U rV in parallel for each spin pair, thus reducing the experi- 
mental intricacies. The diagonal tunneling can be directly im- 
plemented (see Fig lc). The spin-flipping hopping requires 
the even/odd sites to be staggered with A£ st ^10-20kHz, but 
is also efficient (Fig. Id). This scheme, originally developed 
for spin-1 bosons with 3-body interactions [15], leads to the 
interaction picture Hamiltonian in Eq. (1) when the scattering 
is switched off by Feschbach resonances [16]. 

Remarkably enough, starting from this ultracold gas of 
non-relativistic atoms, there are certain regimes where the 
emergent quasiparticles become ultra-relativistic fermions. 
We explore such avenues for a translationally-invariant sys- 
tem, where the bulk energy bands come in degenerate pairs 
£jt ± = £ v 2f v cos £ v cos V ± 2(£ v f 2 sin 2 £ v sin 2 v ) 1//2 , and 
k G [— n, ft] 3 lies in the Brillouin zone. In the 7C-flux regime 
(j) v = n/2, atoms wandering around plaquettes take on an 
overall minus sign, and the bands touch at Nd — 8 differ- 
ent points £ {(d x n , d y n , d z n) : d x ,d y ,d z — 0, l}. Around 
them, low-energy excitations display a relativistic dispersion 
E(p) fa ±(c 2 p 2 + c 2 p 2 + c 2 /? 2 )'/ 2 , where c v = 2t v is the ef- 
fective speed of light, and p = k Aj. Indeed, imposing an 
ultraviolet cutoff \p v \ < 1 /2X C , the effective field theory is 

H^= f dW(r)#*F(r), ff* = £c v a^ v , (3) 

where ^(r) = (ci(r),C2(r),C3(r),C4(r)) f is the field operator, 
CCy = (— \) dv <X v , and p v = —\d/dr v is the momentum. The 
chosen hoppings induce a Clifford algebra {erf , a^} = 2d Vfl , 
and Eq. (3) yields a physical realization of the so-called naive 
Dime fermions in LGT [2]. In our scheme, fermion doubling 
leads to an even number of species which, in contrast to the 
artificial doublers in LGT, correspond to physical and observ- 
able flavors. We remark that each of them has a different chi- 
rality y% = Q5Y5, where 75 = a z g)I, and Q\ = i s 
the axial charge. Chiral symmetry, which plays a fundamental 
role in the standard model classifying right/left-handed parti- 
cles 75 *P = zt'P, cannot be incorporated to the lattice globally. 

Let us stress that we are not limited to the massless limit, 
but can also explore a massive regime with H^ s + H m , where 
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Figure 1: a) Superlattice potential (grey lines) trapping K atoms 
in the main and secondary minima. The hopping between F — 9/2 
levels (in black) is laser-assisted via an intermediate F = 7/2 state 
(in red). The coupling, detuned by d + 8, is induced by an off- 
resonant Raman transition with Rabi frequency fl. b) Scheme of 
the four states of the F = 9/2 manifold (red vertices), connected by 
laser-induced hoppings (green boxes), c) Time-evolution of the pop- 
ulations of the neighbouring hyperfine levels. The solid (dashed) 
line is used for site i (i+l); the red (black) line is used for nif =9/2 
(nif = 7/2). A clear spin-preserving Rabi oscillation between neigh- 
boring sites is shown. The numerical simulation is an exact Runge- 
Kutta time-evolution of the complete model involving all the cou- 
plings and the levels in Fig. la. d) The same as before for a spin- 
ipping hopping. Notice the need for a superlattice staggering (10-20 
kHz) in order to avoid on- site spin-ipping. Exact time-evolution 
shows oscillations between neighboring sites with a different spin. 

H m = J d 3 r v P(r) t mc 2 j3 l I , (r), and J3 = o x ®I. Since no mo- 
mentum transfer is required, this term can be engineered via 
on-site microwave Raman transitions after adiabatic elimina- 
tion of the F = 5 manifold, where mc 2 is the Raman strength. 

In LGT, Wilson envisaged a method to decouple the dou- 
blers from a single Dirac fermion [3]. Our versatile approach 
allows us to realize his idea by combining the previous in- 
gredients with additional tunnelings U rV = — ie~ lipv ' 3 . For 
(p v = n/2, the full Hamiltonian in Eq. (3) turns into 

H? n = Y,c v atp v +m d c 2 l3, m d = m-£(-l) dv m v , (4) 

V V 

where m v c 2 = 2t v depends on the assisted-hopping strength, 
and thus on the laser power. We highlight that the assisted 
hopping strengths scale as f v ~ |i2| 2 /ii. Since |i2| 2 is pro- 
portional to the laser intensities, the tunneling strength, and 
thus the masses, are controlled by the beams power. We fi- 
nally stress that since t v ,t v ~ 0.1-1 kHz, the temperature re- 
quirements of this proposal are of the same order as those 
faced in the context of quantum magnetism in optical lattices, 
currently at the forefront of experimental research. Setting 
mc 2 = 2{t x + t y + t z ), doublers become very massive and de- 
couple from the massless fermion at the center of the Brillouin 
zone Aq = 0. At the expense of breaking chiral symmetry 
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Figure 2: a) In-gap zero-energy modes (dashed red lines) for q = 
(k x ,k y ) = 0, m/4 < m x < 3m/4, and m v = m/2,m z = m/4, for a lat- 
tice with N = 40 3 sites and open boundaries at z = 0,L. b) Boundary 
massless Dirac fermion at z — 0, q = 0, and m x = m y = m z = m/2. 
c) Scheme for a fractional magnetic capacitor consisting of an ax- 
ion well: 6(r,t) = —71 if z £ [zi,Z r ], and 6(r,t) = elsewhere, 
pierced by a magnetic field B = Bz- This is designed by tuning 
m x = m y = m z /2 = m/4 globally, whereas in ^> m is only applied 
to z\ < z < z r - d) Accumulated charge on the "plates" of the capaci- 
tor, for a lattice of N = 30 3 sites, m x = m y = m z /2 = m/4, m = 10m 
(leading to 6 = -% for 12 < z < 18), and flux <p/(f) Q = 2x/15. 



[//p,75] 7^ 0, we have a QS of Wilson fermions invariant un- 
der the anti-unitary operator % a — (il ® O y )<%f, where Jff is 
complex conjugation. 

In an effort to preserve chiral symmetry, domain-wall 
fermions are introduced in 4+1 dimensional LGTs [17], whose 
lower-dimensional descendants are the so-called topological 
insulators [6, 7, 10, 11]. These holographic phases have an 
insulating bulk and metallic boundaries where topologically- 
protected midgap Dirac fermions reside. We can realize these 
phases in experiments by inverting the sign of the Wilson 
mass through the laser intensity. We study the effect of mass 
anisotropy on a lattice with open z-boundaries, which leads 
to the energy spectrum in Fig. 2a. For a critical anisotropy 
m x = 5 m > a mass inversion occurs and some levels leave 
the bulk bands to become zero-energy states; whereas for 
m x = \ m tne midgap states fuse back into the continuum. In 
Fig. 2b, the associated energy band displays a solitary 2+1 
massless Dirac fermion exponentially localized around z = 0, 
whereas its doubler appears at z = L. In other words, there are 
two distant surfaces with an odd number of massless fermions, 
an unambiguous signal of a strong TI that defies the doubling. 

TIs encode a topological order that has dramatic conse- 
quences on their response to electromagnetic fields [6-8]. We 
explore these effects by synthesizing artificial E,B fields [14], 
subjected to a modified Gauss law V E = 4np — (e 2 /nc)V9 ■ 
B, where p is the atomic density. For small masses [6, 7], 
and after a spin rotation, the corresponding Bloch states lead 
to 6 = — jsgn(c x c y c z ) Y*d Qi s g n ( m d)> which only depends on 



the axial charge and the sign of the Wilson mass of each bulk 
Dirac fermion. In Fig. 3a, one observes that the strong TIs 
correspond to a non-vanishing axion field 9 = ^mod(2^). 
Besides, this value is robust with respect to small errors in 
the implementation of the assisted hopping c v = c v + 8c v . 

It is important to notice that the anti-unitary symmetry fixes 
9 = {0, n} [6]. However, it is possible to go beyond this sce- 
nario [18], and examine further possibilities lying in our pro- 
posal, by introducing a complex mass m/3 +im/3ys, leading to 
an arbitrary axion 9 £ [0,27c). Such a term, which breaks 
is obtained via Raman transitions, and in the isotropic limit 
leads to 

9 = -^Le^gnK/) -L^tan- 1 ( — ) , (5) 

L d d \ m dJ 

which contains a perturbation 89 = —Y.d&5 m d tsin ~ l ( m / m d)- 
Note that for m — > 0, 89 — > 0, we recover the previous re- 
sult. Conversely, when m ^S> \mj\, the transparent expression 
89 = —%Y.dQi\ m d\ is obtained (Fig. 3b). This perturbation 
takes on any possible value 89 E [ — 7T, 7r) leading to the to- 
tal axion of Fig. 3c. Therefore, our atomic gas constitutes a 
tunable axion medium where the time- and space-dependence 
of 9(r,t) can be externally adjusted by tailoring the focusing 
width and intensity of the Raman lasers m(r,t). We can now 
explore an exotic consequence of axion electrodynamics: the 
fractional magnetic capacitor. Engineering the axion medium 
in Fig. 2c by focusing a strong Raman laser onto an inner re- 
gion, confines the axion 9 = — 71 to z\ < z < z r , and minimizes 
the effects of the external trapping potential present in experi- 
ments [19]. Besides, we can synthesize a magnetic field gen- 
eralizing the techniques in [20] to a lattice [14]. The field is in- 
troduced via Peierls substitution t v t v e^ ^ drA (likewise for 
f v ) in the Landau gauge A = —(j)yx, where is the magnetic 
flux in units of the flux quantum. Setting (j) = 2np/q, with 
p,q G Z, allows us to diagonalize the open-boundary Hamil- 
tonian in terms of 4gL-bands per momentum inside a reduced 
magnetic Brillouin zone q E [— n, n) x [—n/q,n/q), and then 
obtain the density p(z) from the occupied eigenstates. Then, 
the charge per quantum flux is Q = 2nlgp(z), where lg — 
(c/eB) l l 2 is the magnetic length. Subtracting the value in the 
absence of a magnetic flux leads to Fig. 2d, which confirms the 
continuum prediction Q = Qq8(z — z\) — Qq8 (z — z r ), where 
Qo = | , for an optical lattice with = 30 3 sites, and magnetic 
flux (j) = 2n/15. In such background, the modified Gauss law 
predicts an accumulation of fractional charge per flux quan- 
tum at the boundaries. Therefore, our axion medium plays the 
role of an exotic fractional capacitor whose boundaries act as 
conducting plates. In strong contrast to usual capacitors, the 
charge stored is fractional, and rather than the electric field, 
it is the magnetic field which triggers the effect. Let us recall 
that the 40 K atoms are neutral, and thus this effect corresponds 
to an accumulation of atomic density rather than charge. 

We have described a versatile QS capable of realizing 3+1 
massless/massive Dirac fermions, Wilson fermions, and 3D 
topological insulators as an axion medium. However, we still 
need to address the detection of such effects. To distinguish 
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Figure 3: a) Axion index as a function of the masses m y /m,m x /m, 
and setting m z = m/2. In the invariant regime, only fixed values 
of the axion = {0, n} are allowed, b) Perturbations to the axion 
term 56 in the % -breaking regime, d) Total axion term 6 in the 
^-breaking regime. 

between the different types of relativistic bulk fermions, it 
suffices to measure the number of Dirac points at zero en- 
ergy [21]. This can be accomplished by measuring the atomic 
density as a function of the chemical potential close to zero 
energy [22]; a standard technique that relies on absorption 
images of the Fermi gas (i.e. the shadow projected by the 
atoms on a CCD camera) after the atoms have been released 
from the trap. Considerably more challenging is the detec- 
tion of the edge states characteristic of TIs. The ratio between 
the number of edge and bulk modes makes the direct detec- 
tion by density measurements inefficient. Therefore, new but 
also more demanding methods have been proposed, such as 
the population of edges modes with bosons [19], or the use 
of Raman spectroscopy and Bragg scattering [23] . In this 
Letter, we exploit the consequences of the axion medium to 
propose a density-based measurement. Let us remark that the 
accumulation of density in the fractional magnetic capacitor 
does not suffer from an unbalanced edge-bulk density ratio. In 
fact, an extensive number of atoms will accumulate on the ca- 
pacitor plates, which can be detected by phase-contrast imag- 
ing methods. These methods do not require the trap release, 
such as absorption imaging, but rather recover the atomic den- 
sity in situ by measuring the phase shift of the off-resonant 
light diffracted by the Fermi gas. Accordingly, phase-contrast 
imaging is a non-destructive method that has already been im- 
plemented for bosons and fermions [24]. 

In this Letter, we have presented a feasible scheme of laser- 
assisted tunneling in 3D optical lattices, which allows us to 
design, control, and probe Wilson fermions. This approach 
appears as a promising route towards the first fully-tunable 
realization of 3D TIs. We have shown that fractional mag- 
netic capacitors, predicted by axion electrodynamics, can be 
produced and detected using techniques from optical lattices. 
Besides, since the axion dynamics is tunable and each bound- 
ary can be singled out, phenomena such as the magnetoelec- 
tric effect or the boundary fractional quantum Hall effect can 
also be pursued. Additionally, the interplay between topologi- 
cal order, disorder, and strong interactions can be investigated. 
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